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ABSTRACT 

We discuss two important instability mechanisms that may lead to the limit-cycle 
oscillations of the luminosity of the accretion disks around compact objects: ionization 
instability and radiation-pressure instability. Ionization instability is well established 
as a mechanism of X-ray novae eruptions in black hole binary systems but its ap- 
plicability to AGN is still problematic. Radiation pressure theory has still very weak 
observational background in any of these sources. In the present paper we attempt to 
confront the parameter space of these instabilities with the observational data. At the 
basis of this simple survey of sources properties we argue that the radiation pressure 
instability is likely to be present in several Galactic sources with the Eddington ratios 
above 0.15, and in AGN with the Eddington ratio above 0.025. Our results favor the 
parameterization of the viscosity through the geometrical mean of the radiation and 
gas pressure both in Galactic sources and AGN. More examples of the quasi-regular 
outbursts in the timescales of 100 seconds in Galactic sources, and hundreds of years 
in AGN are needed to formulate firm conclusions. We also show that the disk sizes in 
the X-ray novae are consistent with the ionization instability. This instability may also 
considerably influence the lifetime cycle and overall complexity in the supermassive 
black hole environment. 

Key words: physical processes:accretion; X-rays: binaries; galaxies: active - galaxies: 
evolution - galaxies 



1 INTRODUCTION 

The de scription of the viscous torq ue through the a param- 
eter bv lShakura fc Sunvaev I (|l97of ) boosted the modeling of 
the disk accretion onto the central object, with application 
in various fields, from young stellar objects through stellar 
close binary systems to active galactic nuclei. The param- 
eterization was justified, but not based on a specific well 
developed theory. However, the properties of the stationary 
disk flow weakly depend on this assumption. The spectra 
models were totally unaffected by the viscosity unless either 
departure from the Keplerian motion or departure of the 
emission from a local black body were taken into account. 

However, the stability of disk models does depend sig- 
nificantly on the adopted description of the viscosity. First, 
it was noticed that the assumption of the proportionality 
of the viscous troque to th e total (i.e. gas plus radiation) 
pressure leads to the viscous iPringle et al. (1974) and ther- 
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mal lLightman fc Eardlev I (|l974T ) instabilities. Later, the in- 
stability of the outer gas-pressure dominated parts of the 
disk was discove red, in the region of partially i onize d hydro- 
gen an d helium l|Mever fe Mever-Hofmeister I l|l98lf) , ISmak I 
l| 19841) ). Such instabilities should lead to semi-regular peri- 
odic outbursts. Therefore, observations of objects containing 
accretion disks can be used to test the assumptions about 
the viscosity law. 

Such confrontation of the models and theory is being 
done in case of the ionization instability. The early model 
development was actually motivated by the need to explain 
the dwarf novae outbursts in cataclysmic variables. This in- 
stability is also a leading e xplan ation of the X-ra y novae 



Lasota 



outbu rsts (jCannizzo et al. 1 (| 19821 ): for a rev iew see 

(l200ll) ) . It may also apply to active galaxies dLin fc Shields 
l|l986l) . |Mineshige fc Shields I (QUO)) although observational 
tests are at the early stage of development. 

The presence of the radiation press ure instability is not 
proved yet (see e.g. iDone et al. I ([2007)) . The understand- 
ing of the true nature of the viscosity as being due to the 
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magnetorotational instability iBalbus fc Hawlev I (| 199 if ) did 
not simply set the issue. The limit cycle cannot be seen in 
3-D simulations of this instability since the authors cannot 
follow the global evolution of the disc in a viscous timescale 
and the radial propagation of heating and cooling fronts is 
here neglected. However, recent computations indicate that 
radiation pressure contr ibutes to the viscous torque and the 
instability may be there iHirose et al. I (|2009bh . The compar- 
ison of the predictions of this instability with observational 
data was done only in a few papers so far. 

In this work we propose to study further the models of 
the accretion disc instabilities in a global picture, as well as 
to better support them observationally. In case of the radi- 
ation pressure instability, we use two models: one is based 
on the assumption of the viscoscous torque proportional to 
the total pressure and the other is based on the assumption 
of the viscosous torque proportional to geometrical mean of 
the gas and radiation pressure. We also include the cool- 
ing term due to the outflow. We mark the instability strips 
in the disk radius - acccretion rate plane and we compare 
them with the observed properties of the X-ray lightcurves 
of accreting black holes in binary systems, taken from the 
literature. We claim that the observational constraints for 
the thermal disc instability and the limit-cycle behaviour is 
much more than the one famous example of the microquasar 
GRS 1915+105. We give some examples of other sources and 
discuss further need for the detailed observational studies for 
both Galactic individual X-ray sources and AGNs. 

This article is organized as follows. In Section[2]we dis- 
cuss the theoretical background for the radiation pressure 
and partial hydrogen ionization instabilities. We also present 
the results for the extension and overlapping of the unstable 
regions in the accretion discs. In particular, we focus on the 
constraints for the accretion rate and jet efficiency which 
would be adequate for the astrophysical black hole discs to 
become unstable for one or both types of instabilities. These 
results are based on the numerical codes developed by our- 
selves and discussed in detail in a series of previous works. 
In section [3] we present the observational constraints for the 
disc instabilities found for a number of Galactic black hole 
binaries. We also discuss the supermassive black hole AGNs 
and some observational constraints found in the literature. 
In Section |4] we give a summary and conclusions. 



2 DISC INSTABILITIES 
2.1 Radiation pressure 

The black hole accretion disc with the classical heating 
term proportional to the p ressure with the a coefficient 
|Shakura fc Sunvaev I (|l973l l) is subject to the thermal 
and viscous instability when the radiation pressure dom- 
inates over the gas pressure. This occurs in the inner- 
most radii of the accretion disc around a compact ob- 
ject (in case of a white dwarf such a region cannot be 
present). Radiation pressure instab ility of classical alpha 
models of[S hakura fc Sunvaev I II 1973 ] ') was not i ced ve ry early 
(|Lightman fc Eardlev I (I1974T). iPringle et al. I (|l974l) ) and it 
was fully analyzed bv lShakura fc Sunvaev 1 ll 1976ft . 

The time evolution of the system and its stability is 
governed by the accretion rate outside the unstable region 



(i.e., the mean accretion rate). If the accretion rate is low, 
then the disc remains cold and stable, with a constant low 
luminosity. If the accretion rate is large, then the whole disc 
becomes hot and is stabilized by advecti on, i.e. enters a slim 
disc solution l|Abramowicz et al. I (I1988I D. However, for the 
intermittent accretion rates, larger than some critical value, 
the unstable mode activates. In this case the source enters 
a cycle of bright, hot states, separated by the cold, low lu- 
minosity states. The outburst amplitudes and durations are 
sensitive to black hole mass, viscosity parameter and the 
mean accretion rate. Also, the heating prescription is im- 
portant here. If the viscous heating is proportional to the 
total pressure, the outburst amplitudes are very large, how- 
ever, they can be reduced if the heating is proportional to 
the square root of the gas times the total pressure. If we 
assume the heating proportional only to the gas pressure, 
the instability disappears. In case of the geometrical mean 
of the two components the parameter space of the instability 
is greatly reduced. 

Preliminary shearing-box 3D simulations replacing the 
alpha viscosity with a physical (magnetic) viscosity mech- 
anism indicated that there is no thermal runaway even 
when the ra diation pressure is 1 times larger than the 
gas pressure (IHirose et al. I (12009a! )). However, the same au- 



thors (jHirose et al. ( 2009bl )) in their follow-up work al- 



ready suggested the possibility of radiation pressure insta- 
bility in some of their calculations, as the unstable solu- 
tions may be seen on the surface density - effective tem- 
perature plot. The limit cycle cannot be seen in those 
simulations since they do not follow the global evolution 
of the disc in a viscous timescale and the radial prop- 
agation of heating and cooling fronts is here neglected. 
Full time-dependent computations of the global evolution 
can be only performed with a simple viscosity parameter- 
ization, and in such computations a limit-cycle behaviour 
is seen, wi th disc alternatin g betw e en the hot a n d col d 
states (e.g. iNavakshin et al. I rt200oT). IJaniuk et al. I l|2002l) . 
IJaniuk fc Czerny I (|2007h . ICzernv et al. I l|2009] ))~ 



Observatio nally, the situation is far from clear (see e.g. 
the review by iDone et al. I (j2007)). Several authors sug- 
gested that the radiation pressure instability is an attractive 
explanation of the regular outbursts lasting a few hundreds 
of seconds ob s erved i n the microqua s ar GR S 1915+105 (e.g. 
iTaam et al.~l (jl997l ). iDeegan et ai~l (|2009l )). The radiation 
pressure instability is the only model which explains the ab- 
sence of the direct transitions from the state C to the state 
B in this source. The lightcurves of some other objects also 
show the fluctuations in a form of a limit cycle on the ap- 
propriate timesc ales. An interestin g example is X-ray pulsar 
GRO J 1744-28 l|Cannizzo I (| 19961 ). ICannizzo I (|l998l )) with 
periods of very high accretion rate and low magnetic field 
which allows for the presence of the inner radiation-pressure 
dominated part of the disk. This instability operates when 
the radiation pressure is important, so it is expected only in 
high Eddington ratio Galactic sources and AGN. However, 
for many sources accreting at very high rates the limit cycle 
oscillations have not been reported. Comparison of the Ed- 
dington ratios of the stable sources and those showing fast 
(100 - 1000 s) regular outbursts is a key test of the correct 
viscosity parameterization. 
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2.2 Partial Hydrogen ionization 

Another type of the thermal - viscous instability occurs 
due to the partial ionization of hydrogen and helium. 
The unstable zone is typically present in the outer part 
of the disc, where the effective temperature is of order 
of 6000 K. The ionization instab ility was first fou nd by 
iMever fc Mever-Hofmeister I l|l98ll ) and lSmak 1 1| 19841 ) in the 
disks surrounding white dwarfs and it is at present the 
accepted explanation of the dwarf novae and X-ray novae 
outbursts. The numerical mode l s presented in a num ber of 
works (e.g. iDubus et al. I i|200ll ) I Janiuk et al. I l|2004 ); for a 
review, see lLasota I (|200lf )) studied the one or two (1+1) di- 
mensional models of the global evolution of accretion discs 
under the thermal and viscous instability and confirmed the 
possibility of the limit cycle behaviour. 

The hydrogen ionization instability is thus an exam- 
ple of a firm agreement between the observations and the- 
ory. A recent example is SU UMa-type dwarf nova V344 
Lyr observed by Kepler satellite, for which its outbursts 
and superoutbursts have b een modeled with this instabil- 
ity (|Cannizzo et al. I (|201(J )V 

This instability also likely applies to active galax- 
ies (|Lin fc Shields 1 1| 19861 ). iMineshige fc Shields! l|l99oh al- 
though this is still uncertain. The models of the out- 
bursts were further d evelop ed by ISiemiginowska et al. I 
(|l996h . iJaniuk et af~l (120041), and they were applie d 
to statistics of AGN by ISiemiginowska fc Elvis ll 19971) . 



iMenou fc Quataert I l|200ll) and lHameurv et al. I (|2009l ) ar- 



gued that the amplitude of such outbursts will be small but 
the evaporat ion of the inner dis k enhances the amplitude 
consi derably Janiuk et al J ( |2004l ) and prolongs the qiescent 
state lHameurv et al. ( 20091 ). 

In this case the situation is much similar to the above, 
as the disc cycles between two states. The hot and mostly 
ionized state of a large local accretion rate is intermittent 
with a cold, neutral state of a small local accretion rate. 
Again, the quantitative outcome of the model is governed 
by the assumed external (mean) accretion rate, viscosity and 
mass of the central object. 



2.3 Location of the unstable zones 

We calculated the steady-state models of the accretion 
disc structure, for two exemplary values of the black hole 
mass, characteristic for Galactic sources (10 Mq) and AGN 
(1O 8 M0). The models are based on the vertically averaged 
equations for the energy balance between the viscous heating 
and radiative as well as advective cooling and hydrostatic 
equilibrium. The heating term, governed by the viscosity 
parameter a, is in the radiation pressure dominated region 
assumed proportional to either the total pressure, or to the 
square root of the total times the gas pressure. In the gas 
pressure dominated region, located at larger distances, the 
heating is assumed proportional only to the gas pressure. 
We calculate here the vertical profiles of temperature, den- 
sity and pressure, using the opacity tables that cover the 
temperature range relevant for partial hydrogen ionization, 
including the pr e sence of dust and molecules (see details in 
iRozanska et al. 1 1| 19991 )). 

The basic parameter of each stationary model is the 
global (external) accretion rate, through which we deter- 



mine the total energy flux dissipated in the disc at every 
radius r. Once the effective temperature and surface density 
are determined at every disc radius, we find the stable solu- 
tions, i.e. the accretion rates for which the slope of T — E (or 
M — E) relation is positive, and the unstable solutions, with 
the negative slopes. In other words, the "S-curve" is plotted 
locally at a number of disk radii, and we search for the crit- 
ical rh points at which the curve is bending. These points 
limit the maximum and minimum values of accretion rates 
for which at a given radius the disk will be unstable. In turn, 
we determine the range of radii, for which at a given global 
accretion rate the disc is unstable first due to the radiation 
pressure and then to the ionization instability. 

For the latter, the unstable strip is located at the out- 
skirts of the disk. Obviously, if the instability arises at the 
outer disk, the front will then propagate inwards to much 
smaller radii. However, if the disk size is smaller than the 
inner edge of the unstable strip plotted in Figure [T] no ion- 
ization instability outbursts should take place. Our results, 
based on the detailed vertical structure calculations, are con- 
si stent with the simplified formulae given in the Appendix 
of lLasota I l|200ll) . with respect to the inner boundary of the 
ionisation instability strip. The outer edge we determined 
is somewhat larger, due to a different opacity tables in our 
model which include abso rption on molecul e s, e.g . molecular 
hydrogen, as described in lRozariska et alTI l|l999t ). 

In Figure [1] we show the maps of the disc instabilities 
for the two chosen black hole masses, on the plane radius vs. 
global accretion rate (in dimensionless units). In addition, 
we distinguish the two possible stabilizing mechanisms for 
the radiation pressure instability: the heating prescription 
and the possibility of energy outflow to the jet. The latter 
is parameterized by the following function: 



r/jet = 1 



1 + Am 2 



(1) 



and the jet outflow acts as a source of addit ional cooling 
|Navakshin et al. I (|2000l ). I Janiuk et al. I (|2002h 1. 

The jet outflow reduces the size of the unstable zone for 
large accretion rates, as well as limits the instability to oper- 
ate below some threshold maximum rh. This rate is of course 
sensitive to our adopted parameter for the jet strength. The 
Figure Q] shows the case of a very strong jet, with A = 25 
and in this case the limiting accretion rate is about 20% of 
the Eddington rate. For a 10 times weaker jet the limiting 
accretion rate is about 3 times larger. 

The viscosity parameter, a, only moderately affects the 
results for the radiation pressure instability, apart from the 
timescales of the limit cycle. The Figure [T] presents the 
conservative case of a very small viscosity, a — 0.01, for 
which the extension of the unstable zone is small. If the 
viscosity is larger, the unstable zone size is also increased, 
for instance for a = 0.1 the outer radius of the instabil- 
ity is larger by a factor of ~ 1.25. This also means that 
the minimum accretion rate for which the disc is unsta- 
ble, is slightly smaller. However, in this smallest accretion 
rates the unstable zone is very narrow. For a very narrow 
width of the zone, the instability does not work, because the 
front does no t propagate, and we have only marginally sta- 
ble solutions jS zuszkicwi cz fc Miller I |l997)). For instance, 
lHameurv et al. ( 20091 ) find that the heating and cooling 
fronts of the ionization instability do not propagate strongly 
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Figure 1. The extension of the radiation pressure (solid and dot- 
ted lines) and hydrogen ionization (green; dashed lines) unstable 
zones, depending on mean accretion rates (Eddington units). The 
results are for two heating prescriptions: aPtot (blue; solid lines) 
and ct^PgasPtot (red; dotted lines). The crossed regions mark 
the results for a non-zero fraction of jet power, described by Eq. 
Q with A = 25. The black hole mass is M = 1 X 10 s M Q (bottom) 
and M = 10M® (top). The viscosity is a = 0.01. 



enough in their model to account for the large luminosity os- 
cillations in AGN. 

The outburst cycles caused by the ionization instabil- 
ity are very sensitive to the viscosity parameter. As was 
shown already for dwarf novae, the amplitudes consistent 
with observations can be obtained only in the models with 
non-constant viscosity, i.e. a in the hot state must be larger 
than in the cold state. This may also be the case in the X- 
ray binaries, however in AGN the situation may be different 
(see I Janiuk et al. I (|2004 )). 



Thus the fast radiation pressure instability is expected 
to operate for an average accretion rate higher than a certain 
lower limit, mostly dependent on the adopted viscous scal- 
ing. There is an upper limit as well if the outflow is strongly 
increasing function of the Eddington ratio. The ionization 
instability should operate if the disc is large enough to show 
a partially ionized hydrogen zone for a given average ac- 
cretion rate. Confronting the observational constraints with 
model predictions will in turn allow us to find the constraints 
for the viscosity parameterization and the role of the out- 
flow. In the next sections we make a preliminary step in this 
direction. 



3 OBSERVATIONAL CONSTRAINTS 

The Galactic X-ray binary systems are variable in wide 
range of timescales. First, the transient X-ray sources un- 
dergo their X-ray active states on timescales of years. Sec- 
ond, some sources exhibit X-ray periodic variability on 
timescales of months. Third, some of the most luminous 
sources are variable in timescales of tens-thousands of sec- 
onds. Finally, many of X-ray binaries undergo quasi-periodic 
oscillations. Direct comparison of these data and the mod- 
els is not simple: even the identification of the type of the 
variability with the mechanism is not unique. 

In Table Q] we summarize the properties of the exem- 
plary, best studied X-ray binary sources found in the liter- 
ature, which in our opinion may display radiation pressure 
or ionization instability. We list their characteristic variabil- 
ity timescales and amplitudes, estimated Eddington ratios 
and disc sizes, as well as we indicate a possible instabil- 
ity mechanism responsible for the variability, whenever it 
is in agreement with our computations presented in Sec. [2] 
The maximum disc radius is estimated based on the 60% 
of the Roche lobe size, from the simplified formula given by 
iPaczvnski I (|l97ll ). whenever we had the data for the system 
orbital parameters. 

The estimates given in Table [T] should be treated only 
as indications. The information comes mostly from the lit- 
erature and frequently relies on quantitative description. All 
objects in Table [T] display ionization instability since we se- 
lected objects classified as X-ray novae and there is very 
little uncertainty in the establishment of their nature. Most 
of them show large amplitude outbursts lasting days. The 
ratio of the maximum to minimum flux was estimated from 
the peak and the emission level at the end of the outburst 
so it represents the lower limit - only some of the X-ray 
novae have clear detection in the quiescence. Clearly, more 
careful observational analysis is needed to better study the 
amplitude pattern during the ionization instabilities in these 
sources. 

However, from the comparison of the accretion rates 
and disc sizes with our map presented in Figure [T] one can 
already infer some information about the individual sources. 
All but one among the sources are well within the instabil- 
ity strip, so the instability operates as expected. Only one 
source, SWIFT J1753. 5-0127, is at the inner border of the 
instability strip. However our disc size estimation is based on 
the black hole mass. If this mass is larger than current value 
inferred from the mass function, the unstable strip will be 
broader in its accretion disc. The amplitude of the outbursts 
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Table 1. Sample of the black hole X-ray binary sources. AT is the estimated duration of an outburst, and F max / F m i n is its amplitude. 
Rd/Rs is the estimated disc size in Schwarzschild units. The observations were found in the literature and taken from http:/ /xte.mit.edu/| 



Source 


AT 




M/M E dd 


Rd/Rs 


Instability 


Ref. 


A0620-00 


150 days 


300 


10~ 2 - 3 


4.8 X 10 5 


Ioniz. 


14 


GRS 1915+105 


20-100 yrs 


> 100 


0.25-0.7 


6.3 x 10 5 


Ioniz. 


2 


GRS 1915+105 


100-2000 s 


3-20 


as above 


as above 


Prad 


1,3,19 


GS 1354-64 


~ 30 d 


> 20 


0.1-1.8 


1.8 X 10 5 


Ioniz. 


33 


GS 1354-64 


~ 20 s 


1.5-2 


as above 


as above 


p 

r rad 


4,20,21, 35 


XTE J1550-564 


200 d 


300 


~ 0.15 


1.3 X 10 5 


Ioniz 


36 


XTE J1550-564 


~2000 s 


1.5 


~ 0.15 


as above 


p 

r rad 


5, 22 


GX 339-4 


100-400 days 


75 


< 0.05 


1.6 X 10 5 


Ioniz. 


6, 23, 24 


GRO J0422+32 


200 days 


> 30 


0.002 - 0.02 


4.8 x 10 4 


Ioniz. 


7, 25 


GRO J1655-40 


20-100 days 


16 


5 x 10~ 4 - 0.45 


1.0 x 10 5 


Ioniz. 


8, 26, 27 


GRO J1655-40 


0.1-1000 s 


7.5 


as above 


as above 


Rrad 


8, 32 


4U 1543-47 


50 days 


300 


4.5 X 10~ 4 - 0.04 


9.6 X 10 4 


Ioniz. 


9 


GS 1124-684 


200 days 


24 


~ 10~ 4 - - 1.0 


5.2 x 10 4 


Ioniz. 


6,14 


GS 2023+338 


150 days 


> 100 


0.01 - 1.0 


3.8 x 10 5 


Ioniz. 


10, 29, 30 


GS 2023+338 


60 s ? 


500 


as above 


as above 


Prad 


10 


SWIFT J1753. 5-0127 


150 days 


10 


0.03 


2.0 X 10 4 


Ioniz. 


17, 31 


4U 1630-472 


50-300 days 


60 






Ioniz. 


28 


GRS 1730-312 


6 days 


200 






Ioniz. 


11 


H 1743-322 


60-200 days 


100 






Ioniz. 


12 


GS 2000+251 


200 days 


240 






Ioniz. 


6 


MAXI J1659-152 


20 days 


15 






Ioniz. 




CXOM31 J004253. 1+411422 


> 30 days 


> 300 






Ioniz. 


13 


XTE J1818-245 


100 days 


40 






Ioniz. 


15 


XTE J1650-500 


80 days 


120 






Ioniz. 


16 


XTE J1650-500 


100 s 


24 






p 

r rad 


18 



1 Wu et al. 2010; 2 Deegan et al. 2009; 3 Taam et al. 1997; 4 Revnivtsev et al. 2000; 5 Homan et al. 2001; 6 Tanaka & Shibazaki 1996; 7 
van der Hooft et al. 1999; 8 Harmon et al. 1995; 9 Gliozzi et al. 2010; 10 in't Zand et al. 1992; 11 Trudolyubov et al. 1996; 12 Motta et 
al. 2010; 13 Garcia et al. 2010; 14 Esin et al. (2000); 15 Cadolle Bel et al. 2009; 16 Corbel et al. 2004, 17 Soleri et al. 2008; 18 Tomsick et 
al. 2003; 19 Belloni et al. 2000; 20 Kitamoto et al. 1990; 21 Casares et al. 2004; 22 Sobczak et al. 2000; 23 Hynes et al. 2003; 24 Miller et 
al. 2004; 25 Shrader et 1. (1997); 26 van Paradijs 1996; 27 Kolb et al. 1997; 28 Buxton & Bailyn 2004 29 Zycki et al. 1997; 30 Zycki et al. 
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31 



Zhang et al. 2007; 32 Greiner 1994 



33 Brocksopp et al. 2001; 

1999 



34 Osterbroek et al. 1997; 35 Cui et al. 1999; 36 Sobczak et al. 



in this source is rather small compared to the other X-ray 
novae in Table [T] which would be consistent with the narrow 
instability strip due to too small disk radius. Further studies 
of this exceptional source can provide key tests of the exact 
location of the instability zone. 

We also indicated in Table [T] which sources are promiss- 
ing candidates for the presence of the radiation pressure in- 
stability. In selecting them we paid attention to the possible 
detection of exceptionally low frequency QPO or just re- 
ports of outbursts or variability in timescales longer than 
10 s. Any faster variability than 10 s in unlikely related to 
radiation pressure instability and for such fast QPO there 
are other mechanisms under consideration. 

GRS 1915+105 is the most obvious candidate, with 
its semi-regular outbursts in timescales of 100 - 2000 s 
present in several am ong the brighter characteristic states 
G Belloni et al. ] (|200d )1. Those outbursts were already mod- 
eled by several authors as caused by the radiation pressure 
instability. 

However, fast outbursts are apparently present in sev- 
eral other sources. We show a specific example. In Figure [2] 
we show an exemplary lightcurve of GS 1354-64, which we 
extracted from RXTE data archive. Clearly, a periodicity 
of the ~ 20 s outbursts is visible in the data. The profiles 
with a slow rise and fast decay are characteristic for the 
limit cycle oscillations in the radiation pressure instability. 



For comparison, the lightcurve of Cyg X-l in its hard state 
is plotted in the bottom panel of the Figure [2] We see here 
a very stable X-ray emission and no signatures of the cyclic 
outbursts. Other sources were selected at the basis of their 
description in the literature. The selction is thus not com- 
pletely objective or uniform but may serve as a guide of the 
viability of the approach. 

Having devided the Galactic sources into those which 
possibly show radiation pressure instability and those which 
seem stable we can compare the Eddington ratios within the 
two groups. 

Sources with the Eddington ratio below 0.03 are stable. 
Examples are GRO J0422+32, GX 339-4, as well as Cyg 
X-l. Among the unstable sources, the object XTE J1550- 
564 has the smallest Eddington ratio, 0.15. The instability 
seen i n this source, however, is likely marginal. ICui et all 
( 1999) reported 82 mHz oscillations, whith the frequency 
later increasing to a few Hz. Low f requency oscillatio ns in 
this source were recently studied bv lRao et al. I {2010), and 
they report periods varying between 2 and 10 Hz which is 
far too high for the radiation pressure instability. The inter- 
esting transition hin ting for an instability was reported by 
iHoman et"al~l (|200lh . In the MJD 51,254 observation, when 
the source was still very bright, the luminosity suddenly in- 
creased without a change in the color. Whether indeed this 
single transition hints for the radiation pressure instability 
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or not, the source likley defines the lower limit for the radi- 
ation pressure instability to operate. 

In Table [1] we do not see sources which have very large 
Eddington ratio and are stable against the radiation pressure 
instability. As we mentioned above, GRS 1915+105 is a good 
example of showing outbursts even at Eddington ratio close 
to 1 so it seems we have no upper limit for the Eddington 
ratio in the case of radiation pressure instability. 

Thus, observationally, the radiation pressure instability 
should operate between the Eddington ratio 0.15 up to 1 or 
more. Comparing this with the several theoretical possibili- 
ties plotted in Fig. [1] we can draw certain conclusions. 

First, only the viscosity prescription ctyj P gas Ptot is con- 
sistent with the lower limit for the radiation pressure insta- 
bility, as the unstable region then extends from the Edding- 
ton 0.16 up. The prescription aPtot would allow instability 
to operate at too low luminosity. 

Second, too efficient cooling by the jet is also ruled out. 
The cooling operates similarly in both cases of viscosity pa- 
rameterization and stabilize the disc. For the adopted values 
of the jet efficiency parameter, the disc is stable for Ed- 
dington ratio above 0.22, which is clrearly inconsistent with 
observations. Therefore, the parameter A is Eq[T]of the disk- 
jet coupling must be significantly lower than this exemplary 
value of A — 25. However the jet is by no means excluded 
and still can carry a substantial energy, because in the case 
of equipartition between the disk and jet radiation, for the 
Eddington accretion rate rh — 1 the jet coupling constant 
equal to A — 1 would be enough. 

The parameterization aW P gas Ptot has an additional ad- 
vantage of reducing the outburst amplitude in comparison to 
aPtot. Most of the candidate sources for radiation pressure 
instability show rather low to moderate amplitudes, from 
factor 2 to 20. Only one source - GS 2023+338 - shows huge 
outbursts, wit h the factor of 5 br ightenings in timescales 
of 60 seconds, lin't Zand et al. I l)l992h interpreted this short 
timescale variability as caused by variable absorption. The 
behaviour of this source is exceptional and puzzling. 

When studying the instabilities in the supermassive 
black hole environment, we usually cannot directly observe a 
duty cycle of a one single object, since the black hole masses 
are large and the expected timescales are very long. Instead, 
the statistical studies are useful here and we can find an ev- 
idence for the source episodic activity (e.g. ICzernv et al. I 
(2009)). However, the exceptional o bject is NGC 4395, wi th 
the black hole mass of 3.6 x 10 5 M Q (jPeterson et al. I |2005l )). 
In this source, in principle, we could observe the variabil- 
ity due to radiation pressure instability. As was shown by 
ICzernv et al~l l(2009). the outbursts for the central black hole 
of the mass 1O 7 M0, should last below 100 years, so for a 
mass 30 times smaller, the outbursts should last ~ 3 years! 
No such outbursts are observed. However, this fact is actu- 
ally consistent with our expectations, since the Eddington 
ratio in this source is only 1.2 x 10 -3 . The source is thus 
stable with both aPtot and a^/Pg as Ptot mechanisms and 
provides no useful constraints for the parameterization of 
the viscous torque. 

Significant constraints can be obtained from radio 
galaxies. In case of the accretion discs in radio galaxies, the 
Eddington ratios can be estimated e.g. through the corr e- 
lation with the broad line luminosities (|Dai et al. I |2007l )). 
The FR I and FR II sources in this sample have low Ed- 
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Figure 2. Top panel: an RXTE lightcurve of GS 1354-64, ob- 
served on 05/12/1997. Bottom panel: a lightcurve of Cyg X-l, 
observed by RXTE in the hard state, rescaled to the mean coun- 
trate of GS 1354-64. The time bin is 0.1 s in both data. 



dington ratios, of 0.00975 and 0.0096, for FR I and FR II 
sources, respectively. Observations clearly show that these 
sources are stable against the radiation pressure instability 
since they form very large scale radio structures. In partic- 
ular, the central engine of FR II galaxies must be operating 
in a continuous way for millions of years. Fig. [T] shows that 
their stability is consistent with theory if the heating is given 
by Q.yJ PgasPtot- On the other hand, the FSRQ sources with 
compact radio structures tend to have larger Eddington ra- 
tios. These sources may in fact exhibit episodic activity and 
the small siz e of the structure is in dicating a new episode, as 
proposed bv lCzernv et al.l l|2009l ). Therefore, it seems that 
the assumption of the a-i/ P gas Ptot can accomodate the ob- 
servational constraints both for Galactic sources and AGN. 

A typical value of th e Eddington ratio found in the 
SDSS sample of quasars bv lKellv et all l|2010h is 0.05 with a 
scatter of 0.4 dex. This is also large enough for the episodic 
activity caused by the radiation pressure instability. Possi- 
bly, the selection effect is in fact the reason why we detect 
only the sources in the active state: most of the sources in the 
quiescent state are too dim to be detectable. There is also 
a possibility that active galaxies at high Eddington ratios, 
close to 1, are actually stable due to the stabilizing power of 
jet/outflow. This mechanism seems not to work efficiently 
in Galactic sources but the relative jet power in accreting 
sources rises with the black hole mass: 



log L R = 0.6 log L x + 0.8 log M 



(2) 



as was discussed in the context of th e so called 'fundamen - 
tal plane' of t he bla ck hole activity (I Merloni et al. I l|2003l ). 
iFalcke et al~l (|2004T l). Therefore, this effect in AGN can be 
much stronger than in the Galactic sources. 
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Figure 3. RXTE/ASM lightcurves of two X-ray novae. The bot- 
tom panel shows 4U 1543-475 which is the example of classical 
FRED behaviour. The top panel shows a much more complex be- 
haviour of GRO 1655-40: almost symmetric short outburst, fol- 
lowed by extended phase of activity. This is possibly an analog of 
superoutbursts seen in many CV systems 

4 DISCUSSION 

The theory of accretion disks suggests the presence of two 
instabilities: ionization instability and radiation pressure in- 
stability. In the present paper we made a step towards con- 
fronting the theoretical expectations with the observations 
of Galactic sources and AGN. 

The ionization instability is broadly accepted as an ex- 
planation of the X-ray novae phenomenon and it was com- 
pared to the data for galactic s ources by several authors. 
In this paper we tested whether the size of the disk in the 
X-ray novae systems is consistent with the conditions of the 
ionization instability. The source SWIFT J1753.3-0127 is at 
the border of the instability strip, which may explain the 
low outburst amplitude in this source. All other sources are 
located well within the instability strip supporting this out- 
burst mechanism. 

The details of the outbursts, however, are not well un- 
derstood yet. We note here that the interpretation of the 
time profiles of X-ray novae is somewhat complex. Some of 
the outbursts have well understood FRED profile (i.e. fast 
rise and exponential decay), with the sharp luminosity rise 
due to the the ionization instability and an extended wing 
due to the X-ray irradiation of the outer disc. However, in 
many cases an additional 'superoutburst' follows the first 
outburst (see Figure [3} . The possible interpretation of this 
secondary, extended maximum may be that the accretion 
rate from the co mpanion star increases due to some modu- 
lation effect (see lSmakl (|2010T ) for the analysis of the dwarf 
novae superoutbursts, very much similar to the X-ray novae 
presented here). Therefore the classification and quantita- 
tive analysis of the outburst durations due to the ionization 
instability is not very straightforward. 

Interestingly, the famous microquasar GRS1915+105, 



observed in the very high state since its discovery in 1992, 
may in fact still be in such a 'superoutburst' phase. On the 
other hand, the FRED profile is not always seen possibly 
also due to the lack of irradiation of the disc due to unknown 
reasons. In this case, the ionization instability results in a 
short, symmetric profile. An example of such a source can 
be GRS1730-312, with a timescale of 6 days (see Table [TJ. 

In the case of AGN the applicability of the ioniza- 
tion instability is still under discussion. Fortunately, radio 
observations can give direct insight into the activity his- 
tory. The radio maps of several sources show multiple activ- 
ity periods, mostly in t he for m of double-double structures 
|Schoenmakers et al. I ll2000h. ISaripalli fc Subrahmanvan I 
|2009l) . lMarecki fc Szablewski I j200m . Some of those events 
may be due to mergers and significant change of the jet di- 
rection suggests such a mechanism. However, if the jet axis 
does not change, either a minor merger or ionization insta- 
bility are the likely cause. The timescales can be studied by 
analyzing the ages of the structures. 

Some other sour ces in turn show a decay phase 

l|Marecki fc Swobodal ()2010l )). It is quite likely that the 
microquasar, which has been recently discovered to pos- 
sess the outfl o w hig hly dominated by the kinetic power 
jPakull et al. I (2010)) actually also represents such a fading 
source. The analysis of such a complex behaviour requires 
the radio maps with large dynamical range. Nevertheless, 
multiple radio surveys are under way and more observational 
constraints should be soon available. 

The radiation pressure instability model has been 
proposed firs t to m o del the microqua s ar tim e variability 
|Taam et al. I (|l997l ). iNavakshin et al. I (|2000l ^. In case of 
the regular period i c out bursts of GRS 1915+105 (see e.g. 
iFender fc Belloni I (|2004j M. lasting from ~ 100 to ~ 2000 s, 
(depending o n the source mean lu minosity), this approach 
is successful jjaniuk et al. ] d200l l1. No other quantitative 
mechanism has been put forward to explain the observed 
behavior of this object, and only the limit cycle mechanism 
(likely driven by the radiation pressure instability) explains 
the absence of the direct transitions from its spectral state 
C to the state B. 

Still, the question arises, why the microquasar seems to 
remain an exceptional case where the radiation pressure in- 
stability gives an observational signature. In the present pa- 
per we suggest that several other black hole binaries can also 
be promising objects for the radiation pressure instability. 
All the candidate sources have the Eddington ratio above 
~ 0.15. Such a condition is consistent with theory if the 
viscous torque parameterization as av/ P gas ftot is adopted, 
instead of aP to t ■ Small outburst amplitudes in all candidate 
sources (with one exception) also support a^J P gas -Ptot pre- 
scription. The recent 3-D MHD simulations show the con- 
tribution to the st ress from rad i ation p ressure, but likely 
weaker than cvftot l|Hirose et al. I l|2009bl )). Further numeri- 
cal work and observational constraints aimed at finding the 
proper description of the viscous torque should be treated 
as complementary. 

We argue that the radiation pressure instability also 
likely applies to active galaxies. The derived separations be- 
tween outbursts are on order of 10 6 yrs for a 10 8 Mq black 
hole, while the outburst duration is an order of magnitude 
shorter. 

Again the parameterization of the viscosity through 
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a^J Pgas-Ptot seems to be working better than aPtot- Such 
parameterization is consistent with the lack of instability 
in FR II radio galaxies as their Eddington ratio is below 
0.025 - the lower limit for the instability in active galax- 
ies if the a U P gas -Ptot is adopted. This lower limit here is 
muc h lower than in case of stellar mass objects (see also 
e.g.. ISadowskil (2009)1 due to the direct dependence on the 
black hole mass. This model has been applied recently also 
to explain the apparent young ages of the G iga-Hertz Peak 
Spectrum radiosources (|Czernv et al. I (|2009h ). We speculate 
that in the hot state, the luminous core will power a radio 
jet, while during the cold state the radio activity ceases. 
Scaling the timescale with the black hole mass by a factor 
10 8 gives the outbursts durations of 10 2 — 10 4 yrs, and am- 
plitudes are sensitive to the energy fraction deposited in the 
jet. This gives an additional, model-independent argument 
that the intermittency in quasars on the timescales of hun- 
dreds/thousands of years is likely of a similar origin as in 
the microquasars. 

An important, unstudied aspect is a possible interplay 
between the two instabilities. The location of the unstable 
zone is sensitive to the black hole mass (see Fig. [T|). The 
two zones are located much closer to each other in the case 
of the active galaxies than in the case of Galactic systems. 
In Galactic binaries only one - ionization instability - is fre- 
quently operating. For small accretion rates, which allow for 
systematically lower disc temperatures and the ionization 
instability fits well in the disc size, the radiation pressure 
instability will not develop. If both instabilities are present, 
they are separated in the disc by over 1000 Rs which implies 
that they are well separated in time, acting on timescales 
of seconds and tends of days, respectively, so they can be 
modeled independently, as the radiation pressure instability 
oscillations of a short timescale will be just overimposed on 
the high luminosity state. 

In AGN, the situation is different. The partly ionized 
zone is much closer to the black hole. If the two unstable 
zones are very close to each other, the rate of supply of ma- 
terial to the radiation pressure dominated region may be 
modulated on slightly longer timescales, independently of 
the environment changes in the host galaxy. This poses an 
observational challenge to deconvolve the interplay between 
the two instability timescales, as well as the environmen- 
tal effects. Additional modeling effort should thus be un- 
dertaken to better formulate the theoretical expectations of 
these instabilities. 

Our research has a very preliminary character, and 
further work is clearly needed. On the observational side, 
careful search for excess variability in Galactic sources at 
timescales of a few tens of seconds to a few hundreds of 
seconds should be done. The main difficulty will concern 
inventing a proper mathematical description of this excess, 
because the outbursts due to radiation pressure instability 
are not likely to be strictly periodic and appear clearly in 
periodograms. In the case of the disks around supermassive 
black holes, more constraints should come from detailed ra- 
dio maps of compact sources showing reactivation events 
in short timescales of hundreds - thousands of years. Fur- 
ther development of the models, particularly in the case of 
the radiation pressure instability is also needed, and this 
should be done in two ways. First, if the 3-D MHD compu- 
tations as those done by Hirose et al. (2009) with realistic 



boundary conditions can be prolonged to viscous timescales, 
we would see whether the instability develops a limit cy- 
cle without an ad-hoc parametric description of the viscos- 
ity. This approach is computationally challenging and may 
not happen soon. Second, the sources which in our oppinion 
are promissing candidates for the radiation pressure insta- 
bility do not exhibit regular outbursts which implies that 
strong non-linear/non-local phenomena are important. Sev- 
eral such phenomena can be implemented into the current 
parametric codes. Irradiation is certainly important, and 
actually some models for ionization instability incorporate 
them. The disk irradiation may also be important for the 
radiation pressure instability. The parameter a may depend 
locally on the disk thickness, and in addition the dissipa- 
tion may be coupled to the magnetic pressure with a radius- 
dependent significant time delay (Kluzniak, private commu- 
nication). The theoretical lightcurves can be distorted by a 
stochastic process, such as m agnetic dynamo, operatin g on 
the local dynamical timescale |Maver fc Pringlel (|2006T l ; see 
also I Janiuk fc Czerny I (|2007t ) for the additional discusssion 
of this problem in the context of the magnetically coupled 
hard X-ray corona). That process can evolve e.g. accord- 
ing to the Markov chain model. The timescale and possibly 
also the magnetic cell size are governed by the H/R ratio. 
As a result, the outbursts in the limit cycle can be affected 
by the flickering, as the fluctuations propagate to the inner 
disc. On the other hand, in between the outbursts the fluc- 
tuations are more likely to be smeared out. The resulting 
lightcurve and PDS spectrum will depend on the adopted 
magnitude of the poloidal magnetic field and other param- 
eters. Comparison of such better models with observational 
constraints will help in the future to better understand the 
disk dynamical behaviour. 



5 CONCLUSIONS 

Our preliminary survey of model predictions in confronta- 
tion with observational data suggests that the parametric 
description of accretion disk viscosity through a^J Pgas-Ptot 
is a promissing representation in the radiation pressure dom- 
inated disk part. We selected several Galactic sources as the 
candidates which may show the radiation pressure instabil- 
ity, but further research is clearly needed. The same law 
likely applies to AGN, and the support comes from the sta- 
bility of dwarf Seyfert galaxy NGC 4395 and FR I and FR 
II sources. Ionization instability criterion in Galactic sources 
is consistent with the disk sizes. There are only very limited 
constraints for this instability in AGN. 
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